Competition between symmetry breaking and onset of collapse 
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We analyze the symmetry breaking of matter-wave solitons in a pair of cigar-shaped traps cou- 
pled by tunneling of atoms. The model is based on a system of linearly coupled nonpolynomial 
Schrodinger equations (NPSEs). Unlike the well-known spontaneous-symmetry-breaking (SSB) bi- 
furcation in coupled cubic equations, in the present model the SSB competes with the onset of 
collapse in this system. Stability regions of symmetric and asymmetric solitons, as well as the 
collapse region, are identified in the parameter space of the system. 
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Introduction. It is commonly known that the Gross- 
Pitaevskii equation (GPE) furnishes a very accurate de- 
scription of the Bosc-Einstcin condensation (BEC) in rar- 
efied gases of bosonic atoms [l[ . One important applica- 
tions of the GPE is the prediction of Joscphson oscilla- 
tions and spontaneous symmetry breaking (SSB), alias 
macroscopic quantum self-trapping Q , in double- wellpo- 
tentials (DWPs). Experimentally, the oscillations US] 
and self-trapping Q have been demonstrated in BECs 
with repulsive interactions between atoms, see Ref. @ 
for a review. Asymmetric stationary states trapped in 
the DWP are generated by symmetry-breaking bifurca- 
tions from symmetric or antisymmetric states, in case the 
intrinsic nonlinearity is attractive or repulsive, respec- 
tively 0. In fact, this bifurcation was first predicted (for 
the self-focusing nonlinearity) in the model of dual-core 
nonlinear optical fibers Q ■ 

A natural extension of the DWP is a double-channel 
structure in the two-dimensional (2D) geometry, based on 
two potential wells in one direction (x), which are uni- 
formly extended, as parallel troughs, in the perpendicular 
direction (z) HQ. With the attractive cubic nonlin- 
earity, this structure gives rise to double-hump solitons, 
which arc supported by the DWP in the direction of x, 
while being self-trapped along z, as the ordinary matter- 
wave solitons in the cigar-shaped (single-core) trap fioj ] . 
Similarly shaped gap solitons were predicted in the 
double-channel model with the self-repulsive nonlinearity 
and a periodic potential (an optical lattice, OL) acting 
along z 0, [ll[ . If the nonlinearity is strong enough, an 
obvious symmetric/antisymmetric soliton in the double- 
channel setting with the self-attraction/repulsion may 
bifurcate into an asymmetric mode. This was demon- 
strated both in the full 2D models @, [Tl[ and their ID 
counterparts, which replace the single 2D GPE by a 
pair of ID equations with coordinate x, for wave func- 
tions in the two potential troughs, while the tunnel- 
ing between the troughs in the x direction is approxi- 
mated by a linear coupling between the ID GPEs 
For the attractive nonlinearity, the system of ID GPEs 



is identical to the model of dual-core nonlinear optical 
fibers with anomalous group-velocity dispersion, where 
the symmetry-breaking bifurcation and asymmetric soli- 
tons generated by it have been studied in detail [HI, [l3| . 
The difference between the linearly-coupled models with 
the intrinsic attraction and repulsion is in the character 
of the SSB bifurcation, which is subcritical and super- 
critical in the former and latter cases, respectively. The 
subcritical bifurcation gives rise to a narrow region of 
the bistability, where stable symmetric and asymmetric 
solitons coexist fl2| . 

The DWP may also be extended in two transverse 
directions, forming a pair of parallel pancakes-shaped 
traps. The model developed for this setting is based on a 
pair of linearly-coupled 2D GPEs, which give rise to the 
SSB of 2D solitons and solitary vortices, for cither sign 
of the nonlinearity, provided that the model includes an 
in-plane OL potential (otherwise, all 2D solitons are un- 
stable) Q3. 

The analyses of the SSB in solitons, reported in previ- 
ous works, were dealing with the cubic nonlinearity, ex- 
cept for Ref. pjl , which considered a system of coupled 
equations combining the cubic attractive and quintic re- 
pulsive terms - a setting relevant to optics, rather than to 
BEC. It was demonstrated that the bifurcation diagrams, 
accounting for the SSB of solitons in that system, form 
a closed loop, unlike open diagrams generated by the cu- 
bic nonlinearity. An essential aspect of the description of 
BEC in the low-dimensional setting is that the reduction 
of the GPE from 3D to ID transforms the original cubic 
nonlinearity into a nonpolynomial form [l6l [l7j . In the 
case of the self-attraction, the respective nonpolynomial 
nonlinear Schrodinger equation (NPSE) predicts the oc- 
currence of collapse at a critical value of the density, in 
agreement with the underlying 3D cubic GPE [l6|]. The 
objective of this work is to study a hitherto unexplored 
aspect of the SSB in solitons, viz., the competition be- 
tween the SSB and the onset of the collapse, in the frame 
work of a system of linearly coupled NPSEs (in Rcfs 
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the possibility of the collapse of asymmetric solitons in 
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the 2D dual-core model with the attractive cubic nonlin- 
earity and OL potential was mentioned, but not investi- 
gated, as the SSB happened there at much lower values 
of the density than those necessary for the onset of the 
collapse). Below, we derive the NPSE system, and then 
produce a diagram in its parameter space, which reveals 
regions of stable symmetric and asymmetric solitons and 
a collapse area. 

The model. The starting point of the analysis is the 
scaled 3D GPE for the mean- field wave function, ip, 
which describes the BEC in two parallel identical cigar- 
shaped traps separated by a potential barrier: 
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where z is the longitudinal coordinate, W(z) is the axial 
potential (if any), y is directed perpendicular to the plane 
drawn through axes of the parallel traps, x lies in the 
plane, being orientated perpendicular to the axes, and 
2a is the distance between them. The coefficient 2ng > 
accounts for the attraction between atoms, and fl 2 - for 
the transverse isotropic trapping in each "cigar" . 

Our first objective is to reduce Eq. ((T|) to a system of 
linearly coupled ID equations. To this end, we modify 
the approach developed for the single "cigar" [l6| , adopt- 
ing a superposition of two single-core ansatze: 
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Here /i and f 2 are the ID (axial) wave functions in the 
two cores, and cri. 2 are the respective transverse widths. 

We proceed by substituting ansatz ([2]) into the energy 
functional (Hamiltonian) corresponding to Eq. (fTJ , 
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The underlying assumption, that the distance between 
the "cigars" is essentially larger than the radius of the 
transverse confinement in each of them, implies a 2 Q ^> 
1. Due to this condition, the coupling energy, which 
is produced by the overlap of the two components of 
the wave function in ansatz ([2]), if substituted into 
Hamiltonian ([3]), takes the following form, for g\ = 
(j-2 = cr: E coup i = -K f_™ [fi(z)f%(z) + fl(z)f 2 (z)]dz, 
where k = 2 (aft) 2 exp (— a 2 /cr 2 ) .The main contribu- 
tion to -Ecoupi comes from region x 2 ,y 2 < a 2 around 
the midpoint between the "cigars" . In that region, 
the transverse- confinement radius is determined by the 



ground-state wave function of the 2D harmonic oscillator, 
i.e., a = f^ 1 / 2 , hence the coupling coefficient becomes a 
constant, k = 2 (ail) 2 exp (— a 2 £l) . 

Other terms in Hamiltonian are calculated sepa- 
rately for /i and f 2 . The eventual result is (from now 
on, we set 51 = 1 by means of an obvious rescaling) 

Ua- 2 + a 2 n + 2W(z))\Uz)\ 2 
+ \df n /8z\ 2 - go" 2 \f n (z)\ 4 } dz 

/+oo 
[fi(z)fZ(z) + f*(z)f 2 (z)]dz, (4) 
-oo 

where, as usual [l6j |. the z-dependence of a n is neglected. 
The transverse widths are determined by the variational 
equations, 8E/8a n ~ 0, which leads to the same relations 

as in the single-core model: a 2 = Jl — g |/„| 2 . The sub- 



stitution of this into energy functional (|4]) casts it into 
the final form, 
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The system of coupled NPSEs is derived from Hamilto- 
nian ([S]) as i (f n ) t = 8E/Sf* , i.e., 
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Note that the derivation of the equations for /„ from the 
Hamiltonian in the form of Eq. @ , and subsequent sub- 



stitution of expressions cr 2 = y 1 — g \f n \ , leads to the 
same equations (JB]). Equations © conserve energy ([5]) 
and the total norm (number of atoms in the condensate) , 

N = JIT (l/iWI 2 + l/ 2 Ml") dz = N x + N 2 . 

In the low-density limit, g|/i,2| 2 •C 1, Eqs. © reduce 
to fD GPEs, 

* (fn) t = -(1/2) (/„)„ + W(z)f n - g\f n \ 2 f n - K f 3 - n , 

(7) 

This system, with g > 0, and W(z) = 0, is tantamount 
to the model of the dual-core nonlinear optical fiber with 
the anomalous group- velocity dispersion, where the SSB 
of solitons with f\ (z,t) = f 2 (z,t) was studied in detail 
[12J, ll3( ■ In terms of the double-trap BEC model based 
on system ([7]) with the OL potential, W(z) = ecos(2kx), 
the SSB of regular solitons, for o > 0, and gap solitons, 
for g < 0, was studied in Ref. 0. In the case of g < 0, 
the bifurcation actually breaks the antisymmetry of the 
two-component gap solitons with fi (z, t) — —f 2 (z, t). In 
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FIG. 1: (Color online) . Upper panels: norms Ni and iV2 (solid 
and dashed lines) in the course of the evolution in imaginary 
time. Lower panels: final density profiles of Ui(z) and Ui(z) 
(solid and dashed lines), corresponding to the cases shown 
in the upper panels. The nonlinearity coefficient is g — 0.6, 
while the coupling constant is k — 0, 0.06, and 0.12, in the 
left, central, and right panels, respectively. 



the present work, we focus on the model with the self- 
attraction, g > 0, when the difference of the nonlinearity 
in Eqs. ([7]) from the cubic form is important. As said 
above, the most interesting issue is the competition be- 
tween the SSB and the onset of the collapse, which is 
admitted by the NPSE even in the framework of the ID 
description [l6j . 

Results. To specify the normalization for functions 
fx, 2(2)1 we fix the total norm, N1+N2 = 2, while Nx and 
N2 may vary. Then, the full set of stationary solutions 
to Eqs. ([H]), fi,2(z,t) = e _v *[/i i2 (a;), with real functions 
£^1,2(2) &nd chemical potential fi, is parameterized by two 
coefficients, g and k. We integrated Eqs. ^ in imagi- 
nary time, using a finite-difference Crank-Nicolson code 
and imposing condition N = 2 at each step. Initial wave 
functions were taken as Gaussians with norms Nx = 1.01 
and N 2 = N - Nx = 0.99. 

In Fig. [1] we report typical examples of totally asym- 
metric, strongly asymmetric, and symmetric solitons 
found below the collapse threshold (at g = 0.6), for dif- 
ferent values of the linear coupling k. At k — 0, wave 
function fi(z) with initial norm Nx = 1.01 evolves into a 
single-component bright soliton which absorbs the total 
norm, N = 2, while f 2 (z), with initial norm N2 = 0.99, 
decays to zero. At k = 0.06 (weak linear coupling), the 
wave functions evolve into a stable two-component soli- 
ton with a strongly broken symmetry. Finally, the strong 
linear coupling, with k = 0.12, enforces the evolution of 
the soliton into the symmetric form. 

The collapse of wave function fx (z, t) occurs in simula- 
tions of system ((6]) with K = 0atg>2/3, if the evolution 
commences with initial norm Nx = 1.01. This result 
agrees with known properties of the single-component 
NPSE [lfl, assuming that the total norm, N = 2, goes 
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FIG. 2: The phase diagram of the linearly-coupled system of 
NPSEs in the parameter plane. In regions S and A, the sys- 
tem supports, respectively, stable symmetric [Ui(z) = U2(z)] 
and asymmetric [Ui(z) 7^ ^2(2)] stationary solutions. The 
collapse takes place in region C. The inset is a zoom of the 
area at small values of coupling constant k, which shows the 
border between A and C in detail. 



to fx . A suppression of the collapse takes place at finite 
values of k. The results are summarized in Fig. [5J in 
the form of a diagram in the plane of (<?,k). It features 
three regions: stable symmetric and asymmetric solitons 
in S and A, respectively, and collapsing solutions in C. 
An obvious rescaling argument explains that all solutions 
suffer the collapse at any n for g > 4/3, which is twice the 
above-mentioned critical value, g = 2/3, for the single- 
component system with N = 2. It corresponds to the 
symmetric collapse mode, with equal norms Nx = N2 = 1 
in each trap. In accordance with these features, the bor- 
der between regions A and C in Fig. [2] starts, at K = 0, 
with g = 2/3, while the border between S and C asymp- 
totes to g — 4/3 at large values of k. Notice that in the 
inset of Fig. [5] the (A,C) border seems straight but it is 
actually curved, as shown in the main figure. 

The SSB in solitons is characterized by the asymmetry 
parameter [HEl, © = (Nx - N 2 ) / (Nx + iV 2 ) . The com- 
petition between the symmetry breaking and collapse is 
further illustrated in Fig. [3] by plots of O versus g for 
a relatively weak coupling. In the upper panel we fix 
k = 0.05 and compare the usual model based on linearly 
coupled cubic GPEs, Eqs. (J7J, with the linearly coupled 
NPSEs, Eqs. ©. In the lower panel we plot the NPSE 
curves for three values of k. 

The diagrams of Fig. [3] feature a leap from the sym- 
metric configuration with = to the asymmetric one 
with 7^ 0, cf. a similar situation in the 2D model con- 
sidered in Ref. The transition to asymmetric states 
in the present model always happens by a leap, i.e., the 
symmetry-breaking bifurcation is always subcritical, sim- 
ilar to the situation in the coupled equations with the 
self-attractive cubic nonlinearity [l2| . The SSB may be 
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FIG. 3: (Color online). The asymmetry parameter, O = 
(Ni — Nz) I (Ni + N2), as a function of interaction strength 
g. Upper panel: fixed k = 0.05, as obtained from the coupled 
NPSEs (squares) and ID GPEs, i.e., usual equations with the 
cubic nonlinearity. Lower panel: three values of the param- 
eter k; curves obtained with the coupled NPSEs. The lines 
corresponding to the NPSE system terminate at the collapse 
point. 



considered as a first-order quantum (zero-temperature) 
phase transition, with <d playing the role of the order 
parameter of the transition [18[ • The main difference be- 
tween the 6 (<?) plots generated by Eqs. © and Eqs. 
([7]) is the fact that the former system predicts the col- 
lapse of the asymmetric configuration at g « 0.66, while 
the cubic nonlinearity in the latter system does not give 
rise to any collapse. Lastly, it is seen from Fig. [5] that 
asymmetric solitons do not exist at k > K maJC ~ 0.108. 
Accordingly, with these values of the linear-coupling coef- 
ficient the SSB does not occur, and the collapse happens 
directly in the symmetric mode, at the above-mentioned 
critical value, g = 4/3. 



Conclusion. The objective of this work was to explore 
the SSB in solitons described by the system of linearly 
coupled NPSEs. The difference from the soliton bifur- 
cations in coupled cubic equations is that the SSB com- 
petes with the onset of the collapse. The phase diagram 
in the system's parameter plane was produced, reveal- 
ing regions of stable symmetric and asymmetric solitons, 
and of the collapse as well. The symmetry-breaking bi- 
furcation is subcritical, unless the collapse occurs be- 
fore the SSB. It may be quite interesting to extend this 
model for the description of the SSB in a pair of paral- 
lel pancake-shaped traps, described by linearly coupled 
two-dimensional NPSEs. 
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